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Abstract 

We discuss the phenomenology of the lightest neutralino in models where an 
effective bilinear term in the superpotential parametrizes the explicit breaking of 
R-parity. We consider supergravity scenarios where the lightest super symmetric 
particle (LSP) is the lightest neutralino and which can be explored at LEP2. We 
present a detailed study of the LSP decay properties and general features of the 
corresponding signals expected at LEP2. We also contrast our model with gauge 
mediated supersymmetry breaking. 



1 Introduction 

The search for supersymmetry (SUSY) plays an important role in the experimental pro- 
grams of existing high energy colliders like LEP2, HERA and the Tevatron. It will play 
an even more important role at future colliders like LHC or a linear e + e~ collider. So 
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far most of the effort in searching for supersymmetric signatures has been confined to the 
framework of R-parity-conserving |IJ realizations. Recent data on solar and atmospheric 
neutrinos strongly indicate the need for neutrino conversions [|. Motivated by this 
there has been in the last few years a substantial interest in R-parity violating models 0]. 
The violation of R-parity could arise explicitly as a residual effect of some larger unified 
theory ||, or spontaneously, through nonzero vacuum expectation values (VEV's) for 
scalar neutrinos |6|, In the first case there is a large number of unknown parameters 
characterizing the superpotential of these models, so that for simplicity these effects are 
usually studied assuming in an ad hoc way that only a few dominant terms break R-parity 
explicitly, usually only one. 

We prefer theoretical scenarios which break R-parity only as a result of the properties 
of the vacuum ||. There are two generic cases of spontaneous R-parity breaking models. 
In the first case lepton number is part of the gauge symmetry and there is a new gauge 
boson Z' which gets mass via the Higgs mechanism |J. In this model the lightest SUSY 
particle (LSP) is in general a neutralino which decays, therefore breaking R-parity. The 
LSP decays mostly to visible states such as 

xS - fh, (i) 

where / denotes a charged fermion. These decays are mediated by the Z-boson or by 
the exchange of scalars. In the second class of models there appears a physical massless 
Nambu-Goldstone boson, called majoron. The latter arises in SU{2) (g)C/(l) models where 
the breaking of R-parity occurs spontaneously. In this case the majoron is the LSP, which 
is stable because it is massless (or nearly so). It leads to an additional invisible decay 
mode Xi ~ * v + J i which is R-parity conserving since the majoron has a large R-odd 
singlet sneutrino component [II], |TT| . This decay is absent if lepton number is gauged, as 



the majoron is eaten up by a massive additional Z boson. 

Although models with spontaneous R-parity breaking || [10|, |IT| usually contain addi- 



tional fields not present in the MSSM in order to drive the violation of R-parity (expected 
to lie in the TeV range), they are characterized by much fewer parameters than mod- 
els with explicit breaking of R-parity. Most phenomenological features of these models 
are reproduced by adding three explicit bilinear R-parity breaking terms to the MSSM 



superpotential ||12|| . This renders a systematic way to study R-parity breaking signals 
[ [13| , [14], [L5| and leads to effects that can be large enough to be experimentally observ- 
able, even in the case where neutrino masses are as small as indicated by the simplest 
interpretation of solar and atmospheric neutrino data Q . Moreover, R-parity violating in- 
teractions follow a specific pattern which can be easily characterized. These features have 
been exploited in order to describe the R-parity violating signals expected for chargino 
production at LEP II |16| . 
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Here we consider the phenomenology of the lightest neutralino in the simplest and well 
motivated class of models with an effective explicit R-parity breaking characterized by a 
single bilinear superpotential term ||17|| . Apart from the absence of the majoron-emitting 
x5 decays (which are absent in majoron-less models with spontaneous breaking of Im- 
parity) this bilinear model mimics all features of neutralino decay properties relevant 
for our analysis. For simplicity and for defmiteness we consider supergravity scenarios 
where the lightest supersymmetric particle (LSP) is the lightest neutralino. We present 
a detailed study of the LSP decay properties and general features of the corresponding 
signals expected at LEP2. In the following we denote the minimal SUGRA scenario with 
conserved R-parity by mSUGRA. It is well known that in models with gauge mediated 



supersymmetry breaking (GMSB) the lightest neutralino decays [18, 19], because the 
gravitino is the LSP. We therefore also discuss the possibilities to distinguish between 
GMSB and our R-parity breaking model. 

2 The model 

Here we will adopt a supersymmetric Lagrangian specified by the following superpotential 

W = e ab [h^UjH b 2 + hlQ\D H$ + h%L^RjH" - 0^] + e ab e t L?H b 2 , (2) 

where i,j = 1,2,3 are generation indices, a,b — 1,2 are SU(2) indices, and e is a com- 
pletely anti-symmetric 2x2 matrix, with e X 2 — 1- The symbol "hat" over each letter 
indicates a superfield, with Qi, Lj, Hi, and H 2 being SU{2) doublets with hypercharges 
1/3, —1, — 1, and 1 respectively, and U , D, and R being SU{2) singlets with hypercharges 
— |, |, and 2 respectively. The couplings h v , h D , and h E are 3x3 Yukawa matrices, and 
// and 6i are parameters with units of mass. 

Supersymmetry breaking is parametrized by the standard set of soft supersymmetry 
breaking terms 

v soft = MfQfQ" + Mfu*u j + .\/;f/;; + m^ 2 Z"*Z° + MfR*Rj 

+m 2 Hl H?H? + m 2 H2 H?H2 

- "IM3A3A3 + \M\ 2 X 2 + iM'AiAi + h.c. 

+e ab [^Q a t U 3 H h 2 + A%h%Q h fi 3 H a x + A%h%L b R 3 H« 

-BfjiHfHl + B^LtHl] , (3) 

Note that, in the presence of soft supersymmetry breaking terms the bilinear terms 
£j can not be rotated away, since the rotation, that eliminates it, re-introduces an R- 



parity violating trilinear term, as well as a sneutrino vacuum expectation value [17 . This 



happens even in the case where universal boundary conditions are adopted for the soft 
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breaking terms at the unification scale, since universality will be effectively broken at the 
weak scale due to calculable renormalization effects. For definiteness and simplicity we 
will adopt this assumption throughout this paper. 

Although for the discussion of flavour-changing processes, such as neutrino oscilla- 
tions involving all three generations, it is important to consider the full three-generation 
structure of the model, for the following discussion of neutralino decay properties it will 
suffice to assume R-parity Violation (RPV) only in the third generation, as a first ap- 
proximation. In this case we will omit the labels i,j in the superptotential and the soft 



breaking terms [17, 20 



W 

V S0 ft 



h t QsU 3 H 2 + hQ^H, + KL^H, - fiH,^ + e 3 L 3 H 2 



A t h t Q a 3 U 3 H b 2 + A b h b Q\D 3 Hl + A T h T L b 3 R 3 H? 



—BfxH^H 2 + B 3 e 3 L 3 H 2 



+ mass terms. 



(4) 
(5) 



This amounts to neglecting the Jjt p effects in the two first families. 

The bilinear terms in Eqs. (f|) and (|5[) lead to a mixing between the charginos and the 
r-lepton which is described by the mass matrix 



Mr 



M 2 -^9^2 



72^3 



-e 3 





-75M3 
73 Ml 



(6) 



where t>i, v 2 , and v 3 are the vevs of H®, H®, and z/ r , respectively. As in the MSSM, the 
chargino mass matrix is diagonalized by two rotation matrices U and V 

n-± 
m-± 
m T 

The lightest eigenstate of this mass matrix must be the tau lepton (r^) and so the mass 



U*M C V 



-1 



(7) 



is constrained to be 1.7771 GeV. As explained in ||21|| , the tau Yukawa coupling becomes 
a function of the SUSY parameters appearing in the mass matrix. 
The neutralino mass matrix is given by: 



(8) 



where gi = g'/2 and g 2 = g/2 denote gauge couplings. This matrix is diagonalised by a 
5x5 unitary matrix N, 





" M l 





-givi 


91V2 


-9iv 3 ~ 







M 2 


92V1 


-92V2 


92V3 


M N = 


-QiVi 


92V1 












91V2 


-g2v 2 







£3 




. -91V3 


92V 3 





£3 






xi 







(9) 
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where V° = (~iB, -iW 3 , H d , H u , u T ). 

The up squark mass matrix is given by 



Mi 



Mq + \v 2 h 2 + A UL 



(v 2 A u - fl Vl 
and the down squark mass matrix by 



M 2 + \v\h 2 + A 



UR 



d 



Mq + \v\h d 2 + A DL 



M|, + \v 2 h d 2 + A DR 



(10) 



(11) 



— l ( n 

dl - g{- g 

y\v 2 -v 2 + v 2 ). 



W 2 )(v!-v 2 + vi), A UR 



The sum of the v 2 is given by 



with A UL = \{g 2 
\g' 2 (vl -v\ + «f), and A DR 
m w = 9 2 ( v i + v 2 + v l)/2- The mass eigenstates are given by §1 = cos #~+ q R sin #~ and 
92 = QrcosO- — qL sin 0~. The sfermion mixing angle is given by 



cos 6~ 



-M 2 

Qlr 



sin 8~ 



M 2 

ILL 



91 



(M 2 -m 2 ) 2 + (M 2 f q J(M 2 -m 2 ) 2 + (M 2 



11' 



Qlr' 



(12) 



111 



Qi' 



Qlr' 



Here M~ are the corresponding entries of the mass matrices in Eqs. (0) and (|iT|). 

In addition the charged Higgs bosons mix with charged sleptons and the real (imag- 
inary) parts of the sneutrino mix the scalar (pseudoscalar) Higgs bosons. The formulas 

I and are reproduced, for completeness, in the appendix. The cor- 



can be found in 21 



responding mass eigenstates are denoted by for the charged scalars, S® for the neutral 
scalars, and P° for the pseudoscalars. 



3 Numerical results 

In this section we present numerical predictions for the lightest and second lightest neu- 
tralino production cross sections in e + e~ collisions, namely, e + e~ — > XiXi^ XiX2- Moreover 
we will characterize in detail all branching ratios for the lightest neutralino decays, which 
violate R-parity. 

The relevant parameters include the ft p parameters and the standard mSUGRA pa- 
rameters M1/2, m o, tan/?, where M1/2 is the common gaugino mass, m the common 
scalar mass, and tan/? = Vijv\ is the ratio of the vacuum expectation values of the Higgs 
fields. The absolute value of /1 is fixed by radiative breaking of electroweak symmetry. We 
take n positive to be in agreement with the b — > 57 decay |p3 |. As representative values 



of tan/3 we take tan/5 = 3 and 50. It is a feature of models with purely spontaneous 
breaking of R-parity that neutrinos acquire a mass only due to the violation of R-parity 
|], |7|, This feature also applies to models characterized by purely bilinear breaking of 
R-parity, like our reference model charaterized by Eqs. (Q) and (^). As a result the ft p 
violating parameters are directly related with m U3 , the mass of the neutrino u 3 , which is 
generated due to the mixing implicit in Eq. (ffl). 
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3.1 Neutralino Production 



While the violation of R-parity would allow for the single production of supersymmetric 
particles [jn|, for the assumed values of the Jft p violation parameters indicated by the 
simplest interpretation of solar and atmospheric neutrino data @, |3| , these cross sections 
are typically too small to be observable. As a result neutralino production at LEP2 in 
our model typically occurs in pairs with essentially the same cross sections as in the 
mSUGRA case. In Fig. [I|a and b we show the maximum and minimum attainable values 
for the e + e~ — > XiXi an d e + e~ — > X1X2 production cross sections as a function of m^o at 
\/s = 205 GeV. We compare the cases tan/? = 3 and tan/? = 50, varying M1/2 between 
90 GeV and 260 GeV and mo between 50 GeV and 500 GeV. One can see that, indeed, 
these results are identical to those obtained in the mSUGRA. The X1X1 production cross 
section can reach approximately 1 pb. In our calculation we have used the formula as given 
in |25] and, in addition, we have included initial state radiation (ISR) using the formula 



given in [ 26[] . Note that e l and e# are exchanged in the t- and w-channel implying that a 



large fraction of the neutralinos will be produced in the forward and backward directions. 

In order to show more explicitly the dependence of the cross sections on the parameters 
m and M1/2 we plot in Fig. |^a and b the contour lines of <r(e + e~ — > X1X1) i R the m®- 
M1/2 plane at y/s = 205 GeV for tan (3 = 3 and tan/? = 50. The contour lines for 
a(e + e _ — > X1X2) are given in Fig. ^|c and d. 



3.2 Neutralino Decay Length 

If unprotected by the ad hoc assumption of R-parity conservation the LSP will decay 
as a result of gauge boson, squark, slepton and Higgs boson exchanges. The relevant 
contributions to these decays are given in Table 1. The Feynman diagrams for the decays 
not involving taus, i.e. Xi ^3// (/ = e > p ei A 4 ; ^> u , d, c, s, b) are shown explicitly in 
Fig. |. 

For the observability of the R-parity violating effects it is crucial that with this choice 
of parameters the LSP will decay most of the time inside the detector. The neutralino 
decay path expected at LEP2 depends crucially on the values of ]fl p violating parameters 
or, equivalently, on the value of the heaviest neutrino mass, m U3 . We fix the value of 
m Vi as indicated by the analyses of the atmospheric neutrino data ||. It is important to 
note that, as explained in due to the projective nature of the neutrino mass matrix 
[0], only one of the three neutrinos picks up a mass in tree approximation. This means 
that, neglecting radiative corrections which give small masses to the first two neutrinos 
in order to account for the solar neutrino data, the neutralino decay length scale is set 
mainly by the tree-level value of m U3 . In ref. we have explicitly shown that this is a 
good approximation for most points in parameter space. 
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Figure 1: Maximum and minimum attainable values for the e+e — > X1X1 (full lines) 
and e+e - — > X1X2 (dashed lines) production cross sections in fb as a function of m^o for 
y/s = 205 GeV, 50 GeV < m < 500 GeV, 90 GeV < M 1/2 < 270 GeV, a) tan/3 = 3, and 
b) tan (3 = 50. ISR corrections are included. 

In Fig. H we plot the Xi decay length in cm expected at LEP2 for y/s = 205 GeV. Here 
and later on we consider the neutralinos stemming from the process e+e~ — > X1X1 when 
discussing the decay length. In Fig. |]a we plot the Xi decay length in cm as a function 
of neutrino mass m V3 , for different m^o between 60 and 90 GeV, with m = 100 GeV, 
and tan/3 = 3. As can be seen the expected neutralino decay length is typically such that 
the decays occur inside the detector, leading to a drastic modification of the mSUGRA 
signals. An equivalent way of presenting the neutralino decay path at LEP2 is displayed 
in Fig. ^b, which gives the decay length of Xi as a function of m^o for m U3 = 0.01, 
0.1, and 1 eV. Finally, we show the dependence of the neutralino decay path on the 
supergravity parameters fixing the magnitude of Tjl v violating parameters or, equivalently, 
the magnitude of the heaviest neutrino mass, m U3 . In Fig. [5]a and b we plot the contour 
lines of the decay length of Xi i n the m Q -Mi/ 2 plane for m Vi = 0.06 eV, tan/3 = 3 and 50. 
Note that the decay length is short enough that it may happen inside typical high energy 
collider detectors even for the small neutrino mass values ~ 0.06 eV indicated by the 
atmospheric neutrino data 0. For large values of tan/3 the total decay width increases 
and, correspondingly, the decay path decreases due to the tau Yukawa coupling and the 
bottom Yukawa coupling. 
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Figure 2: Contour lines of the production cross sections in fb, in the m -M 1 / 2 plane for 
= 205 GeV, a) e+ e - - x°iX°i, tan/3 = 3, b) e + e~ - x?Xi, tan/3 = 50, c) e+e~ - x\xl 
tan/3 = 3, and d) e + e~ — > tan/3 = 50. ISR corrections are included. 
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Figure 3: Feynman graphs for the decay Xi ~^ ^3 f f where / 7^ r. 
3.3 Neutralino Branching Ratios 

As discussed in the beginning of this section, the lightest neutralino Xi will typically decay 
in the detector. In the following we present our results for the branching ratios of all Im- 
parity violating 3-body decay of x?, and of the radiative decay Xi ~^ v 3l- The Feynman 
diagrams for the decays Xi v zf f (f — e > u e, fJ>, Vp, u, d, c, s, b) are shown in Fig. [3[ For 
this class of decays we have Z°, Pf, and S® exchange in the direct channel (Fig. ^a and b) 
and / exchange in the crossed channels (Fig. [3]c and d). In particular in the case / = b the 
Pj° and Sj exchange contributions are significant. This is quite analogous to the results 
found in |27| for X2 ~^ Xiff decays. The particles exchanged in the s-, t-, and ^-channel 



for the decays Xi ~^ r± l T V[, (I = e,/i), Xi ~~ * ^qq' (q,q f = u,d,s,c), Xi ~^ i~~t + Vi, and 
Xi — ► 3^3 are given in Table [I] . 

In the calculations we have included all mixing effects, in particular the standard 
MSSM /x — fa mixing effects and those induced by the bilinear R-parity violating terms, 
i.e. Re{y T ) - h° - H°, Im(h) - A - G°, [[§], R - H± - G ± m, u T - xl and r 



Xj mixings [T3|. These mixing effects are particularly important in the calculations of 



the various R-parity violating decay rates of Xi, which are discussed below. 

In the following plots Fig. |6| - |13] we show contour lines in the m^-Min plane for 
the branching ratios in % of the various Xi decays, in (a) for tan/3 = 3 and in (b) for 
tan/? = 50. We have fixed the mass of the heaviest neutrino to m V3 = 0.06 eV 0. It 
turns out, that in the range 10~ 2 eV < m V3 < 1 keV all the X® decay branching ratios 
are rather insensitive to the actual value of m U3 . This is an important feature of our 
supergravity-type R-parity violating model. It is a consequence of the fact that, as a 
result of the universal supergravity boundary conditions on the soft breaking terms, all 
R-parity violating couplings are proportional to a unique common parameter which may 
be taken as es/fi. For a more detailed discussion on this proportionality the reader is 
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Figure 4: Decay length of the lightest neutralino in cm for ^/s = 205 GeV, a) as a function 
of m U3 for m^o = 60, 70, 80, and 90 GeV, b) as a function of m^o for m U3 = 0.01, 0.1, and 
1 eV. 
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Figure 5: Decay length of the lightest neutralino in cm in the m -Mi/ 2 plane for y/s = 
205 GeV, a) tan (3 = 3, and b) tan/3 = 50. The R-parity violating parameters are fixed 
such that m„„ = 0.06 GeV. 
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Figure 6: Branching ratios for Xi ~ y 3 v in % in the rriQ—Mi/2 plane for a) tan/5 = 3, and 
b) tan/? = 50. The R-parity violating parameters are fixed such that m U3 = 0.06 GeV. 

referred to ref. ||. Also note that for M1/2 > 220 GeV the neutralino mass becomes 
larger than m w and so that xi decays into real W and Z are possible. The effects of 
these real decays can be seen for M1/2 > 220 GeV in most of the following plots. For the 
large tan (3 case (tan (3 = 50) and M1/2 3> M the mass of the lighter charged boson S 1 is 
smaller than m^o (upper left corner of Fig. ^b - p~3|b) . In this region of the parameter space 
the two 2-body decays Xi —> W ± t± and Xi ~ * Sf T ± compete. The first one is R-parity 
violating, but has more phase space than the second one which is R-parity conserving, 
since is mainly a stau. For this reason, the most import final state is t + t~z/3, followed 
by r ± gg / and r^l^Ui (I = e,/i) as shown in Figs. [12|, [10], and §, respectively. All other 
final states have nearly vanishing branching ratios in this corner of the parameter space. 

Fig. |B|a and b exhibit the contour lines for the branching ratio of the invisible decay 
Xi —>■ This branching ratio can reach 7% for the parameters chosen. In Fig. |7| 

and [8] we show the branching ratio for the decays Xi ~ > u 3 l + l~ an d Xi ~ y v^qq where 
/ and q denote the leptons and quarks of the first two generations, summed over all 
flavors. These branching ratios can go up to 3% and 15%, respectively. Notice that the 
sneutrino, slepton, and squark exchange contributions to the Xi decays become larger 
with increasing mo, despite the fact that the increase of the scalar masses m^, mj, trig 
suppresses these exchange contributions. This trend can also be observed in Fig. §, || 



and 10. This happens because the tadpole equations correlate /i to m . Increasing fx 
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Figure 7: Branching ratios for Xi ^ ^3 / + / in % in the rriQ— M\m plane for a) tan/? = 3, 
and b) tan/3 = 50. Here I is the sum of e and fi. The R-parity violating parameters are 



fixed such that m, 



0.06 GeV. 



while keeping M\ and M 2 fixed implies increasing the gaugino content of Xi an d, hence, 
enhancing the Xrf-f couplings. 



In Fig. |9] and |10| we show the contour lines for the branching ratios of the LSP decays 
involving a single tau, namely Xi ~^ vit^P and Xi ~ * T± 10* '1 where I, q, and q' are 
summed over the first two generations. The branching for these decay modes can reach 
up to 20% and 60% respectively. For Mi/ 2 > 220 GeV decays into real W ± dominate. If 
this is the case and if both x\ produced in e + e~ — > X1X1 decay according to these modes 
this would lead to very distinctive final states, such as 4jr + r + , r + r + / - / - (/ = e,/z), 
or T + r + e~fi~. The full list of expected signals is given in Table |2] . The first column 
in this table specifies the two pairs of Xi decay modes, while the second one gives the 
corresponding signature. In the last column we state whether the corresponding signature 
exists for e + e~ — > X1X2 production within mSUGRA. 

The LSP decays involving only third generation fermions, namely, Xi Vsbb and 
Xi — > v-z t + t~ are different from those into the first and second generation fermion pairs, 
because the Higgs boson exchanges and the Yukawa terms play a very important role. 



This can be seen in Fig. [11] and [12] , where we plot the contour lines for these decays. The 
branching ratio of Xi ~ * Vzbb can reach up to 97%. The decay rate is large because the 
scalar exchange contributions (S®, P®, b k ) are large for Mi/ 2 < 200 GeV. Note that this is 
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Table 1: Contributions involved in the lightest neutralino 3-body decay modes. The s-, 
and it-channels are defined by: s — {p\ — P2) 2 , t — (pi — P3) 2 , and u — (pi — P4) 2 . See 
also Fig. [| 
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Figure 8: Branching ratios for — > u 3 q q in % in the m^-Mx^ plane for a) tan j3 = 3, and 
b) tan (3 = 50. Here q is the sum over u, d, s, and c. The R-parity violating parameters 



are fixed such that m 



0.06 GeV. 



also the case for tan (3 = 3, because not only the neutrino-neutralino mixing proportional 
to m U3 is important but also the neutrino-higgsino mixing proportional to e^/fi. The 
decrease of the branching ratio with increasing m is due to the decrease of the higgsino 
component of Xi an d the increase of the Higgs boson masses. For M1/2 > 200 GeV the 
decays into real W + and Z° are possible, reducing the branching ratio of Xi ~^ v$bb. 
As shown in Fig. O the branching ratio for Xi ~> ^3 T+ T ~ is very small for tan j3 = 3 
and M1/2 200 GeV. This is due to the destructive interference between Z° contribution 
and the contributions of the exchanged charged scalar particles (mainly due to the stau 
components of S^). 



Finally we have also considered the radiative LSP decay mode Xi 



^3 7 



In 



Fig. [13] the branching ratio for this mode is shown. This decay proceeds only at one-loop 
level and therefore is in general suppressed compared to the three-body decay modes. 
However, for My 2 < 125 GeV and large tan/? it exceeds 1%, leading to interesting sig- 
natures like e + e~ — > X1X1 ~ > T± ^ T 1 + $t- Due to initial state radiation it can easily 
happen that a second photon is observed in the same event. 

The complete list of possible signatures stemming from LSP decays in our bilinear Jji v 
model is shown in Table ||] . In this table we also indicate whether the same signatures 
could also arise in mSUGRA as a result of e + e" -> X1X2 followed by the MSSM decay 
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Figure 9: Branching ratios for Xi ~ > ^ r± ^ T i n % in the rriQ-Mi/2 plane for a) tan/? = 3, 
and b) tan/3 = 50. Here I is the sum of e and fi. The R-parity violating parameters are 



fixed such that m h 



0.06 GeV. 



modes of x° if its production is kinematically allowed. The final states 4 jets + r + 
2 jets + $t, and r + (e or n) + ^ r would also occur in mSUGRA via the decay of x° 
into xf. However one expects in general that these decay modes are suppressed within 
mSUGRA. In contrast in the R-parity violating case these signatures can be rather large 
as can be seen from Fig. [5] and [K]. Note moreover, that some of the signatures are 
practically background free. For example, due to the Majorana nature of x\i one can have 
two same-sign r leptons + 4 jets + ^/t- Other interesting signals are: r + 3 (e and/or 
/i) + tf T , 3 r + (e or fi) + j/ T , r + (e or yu) + 2 jets + tf T , r + 4 jets + rf T , t^t^ + (e 
or /i) + 2 jets + ^t, or r ± r ± + / T / /=F + ^ with / = e, /i. In Table [| we give masses and 
branching ratios for typical examples. 

As it is well known, also in gauge mediated supersymmetry breaking models (GMSB) 
Il| the neutralino can decay inside the detector, because the gravitino G is the LSP. It 



is therefore an interesting question if the R-parity violating model can be confused with 
GMSB. To answer this question let us have a look at the dominant decay modes of the 
lightest neutralino in GMSB. If the lightest neutralino is the NLSP, its main decay mode 
in GMSB is 
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Figure 11: Branching ratios for Xi ~* v z bb in % in the mo-Mi/ 2 plane for a) tan/3 = 3, 
and b) tan (3 = 50. The R-parity violating parameters are fixed such that m U3 = 0.06 GeV. 
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Combination of \i decay modes 


signature 


mSUGRA-like 


(3 1/) (3 i/) 




yes 


(3 i/) (z/ 3 /+/-) 


2 leptons + fa 


yes 


(3 1/) (u 3 qq) 


2 jets + ^ T 


yes 


(3 i/) (i/ 3 66) 






(3 1/) (z/^r ± / T ) with I — e, fx 


r + (e or /i) + 


yes, but suppressed 


(3 1/) (r^g') 


r + 2 jets + fa 


yes, but suppressed 


(3z/) (1/37) 


1 + fa 


yes 


(z/ 3 Z+Z-) (*/ 3 z /+ r) 


4 leptons + fa 


no 


(i/ 3 Z+Z _ ) (z/ 3 gg) 


2 leptons + 2jets + fa 


no 


(z/ 3 /+/-) (1/366) 






{v 3 l + l~) (uitH^) with / = e,/x 


t + 3 (e and/or //) + ^ r 


no 


(z/ 3 t + t~) (^r ± / T ) with I = e, fx 


3 r + (e or fx) + fa 


no 


(i/ 3 Z+Z ) (r ± g? / ) 


r + 2 leptons + 2 jets + $r 


no 


(i/ 3 r+r-) (r^') 


3 r + 2 jets + fa 


no 


(i/ 3 Z+Z-) (1/37) 


2 leptons + 7 + fa 


no 








(i/ 3 gg) (^366) 


4 jets + fa 


yes, but suppressed 


(1/366) (1/366) 






(z/ 3 gg) (vitH*) with I — e, fx 


r + (e or fx) + 2 jets + ^ T 


no 


/ 7 7 \ / 1 * —1— \ • 7 7 

(1/360) [viT l + ) with I — e, fx 






(u 3 qq) (T k qq') 


r + 4 jets + $r 


no 


(1/366) (r^qq') 






(vsqq) (^37) 


2 jets + 7 + fa 


no 


(1/366) (1/37) 






far^) (^r±Z' T ) 


T^T* + Pl' T + fa 


no 




f^T* + Z^Z'* + fa 


no 


(i/ i r ± ZT) (r±^) 


r r + (e or /x) + 2 jets + $r 


no 




r ± r T + (e or yu) + 2 jets + 


no 


(^Z*) (i/37) 


r + (e or /i) + 7 + 


no 


(r ± gg / ) (r^g') 


+ 4 jets + $t 


no 




+ 4 jets + $r 


no 


(r ± gg / ) (i/37) 


r + 2 jets + 7 + $r 


no 


(^37) (^37) 


2 7 + fa 


no 



Table 2: The signatures expected from the process e + e — > X1X1 m the bilinear x% p model. 
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0.6 


5.6 


1 8 


0.5 


1.8 


17.8 


BR(x? - gg'r±) 


1.1 


16.1 


53.7 


0.9 


5.1 


53.2 


BR(x? - 66z/ 3 ) 


96.5 


62.6 


13.4 


97.1 


88.4 


13.3 


BR(x? - r-r+z/ 3 ) 


0.1 


1.5 


8.6 


0.5 


1.0 


8.8 



Table 3: Masses and branching ratios for the points: A (Mi/2,?tiq) = (153,155), B 
(Mi/2, mo) = (153,440), and C (M 1/2 ,m ) = (251,440) for both tan/5 = 3 and 50. The 
masses are given in GeV and the branching ratios in % and we only give those larger than 
0.1%. Here the same summations of the final states are performed as in the figures, nig 
is the averaged squark mass for the first two generations. 

where G is the gravitino. For the case where at least one of the sleptons is lighter than 
the lightest neutralino the latter has the following decay chain Xi ~> I ^ T —>■ l T G ,. 
In principle three-body decay modes mediated by virtual photon, virtual Z-boson and 
virtual sfermions also exist. However, in the neutralino mass range considered here these 
decays are phase-space-supressed |L8|, |3(J. This implies that the R-parity violating model 



can not be confused with GMSB, because (i) in GMSB the final states containing quarks 
are strongly suppressed, and (ii) GMSB with conserved R-parity implies lepton flavour 
conservation, and therefore there are no final states like e + e + r~T~ + $t- A further 
interesting question would be how the neutralino phenomenology changes in a GMSB 
scenario with broken R-parity. The main consequence would be an enhancement of final 
states containing photons and/or leptons. A detailed study of this question is, however, 
beyond the scope of the present paper. 
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4 Conclusions 



We have studied the production of the lightest neutralino Xi a ^ LEP2 and the result- 
ing phenomenology in models where an effective bilinear term in the superpotential 
parametrizes the explicit breaking of R-parity. We have considered supergravity sce- 
narios which can be explored at LEP2 in which the lightest neutralino is also the lightest 
supersymmetric particle. We have presented a detailed study of the LSP Xi decay prop- 
erties and studied the general features of the corresponding signals expected at LEP2. 
A detailed investigation of the possible detectability of the signals discussed in Table || 
taking into account realistic detector features is beyond the scope of this paper. Clearly, 
existing LEP2 data are already probing the part of the parameter region which corre- 
sponds to approximately m^o < 40 GeV. Finally, we note that, in addition to important 
modifications in the x? decay properties, R-parity violating decay models lead also to new 
interesting features in other decays, such as charged plj and neutral |28| Higgs boson and 
slepton decays, stop decays |31], |22|, and gluino cascade decays In addition we 



have shown that the R-parity violating model can not be confused with gauge mediated 
supersymmetry breaking and conserved R-parity due to the absence of several final states 
in the GMSB case. 
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A Scalar Mass Matrices 



The mass matrix of the charged scalar sector follows from the quadratic terms in the 
scalar potential ET|, |22 . 

y quadratic = S' M|±S' + (13) 

whereS' - = [H^ , H~2 , f£ , Tr]. For convenience reasons we will divide this 4x4 matrix 
into 2x2 blocks in the following way: 



M 2 



s± 



2 T 



(14) 
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where the charged Higgs block is 



(15) 



Bn + \g 2 viv 2 



and h T is the tau Yukawa coupling. 



M? 



ti 



By + z g 2 v iv 2 
By%+ 1 1 g 2 (v 2 + v 2 )-B 3 e 3 % + % 



' V2 



(< 2 



(16) 



\h\v\-\g\v\ 



^ a _ e a + a 



m% 3 + \h%(v\ + vl) - \g' 2 (v\ - v 2 + v 2 ) 



The mixing between the charged Higgs sector and the stau sector is given by the following 
2x2 block: 



Mb* 



-fjLe 3 - \h 2 T v x v 3 + \g 2 v x v 3 -B 3 e 3 + \g 2 v 2 v 3 
-^K(e 3 v 2 + A T v 3 ) —^=hr(fjLV3 + e 3 v ^ 



(17) 



As we see the charged Higgs bosons mix with charged sleptons. 

In a similar way the real (imaginary) parts of the sneutrino mix the scalar (pseu- 
doscalar) Higgs bosons. The quadratic scalar potential responsible for the neutral Higgs 
sector mass matrices includes 



±P'° T M?, P' + S' 0X MLS'° + . . . 



/0T T 



Vquadratic 2 1 ivl po 1 ••• o lu s ' 
)/0 T r „0 ~Il o/0 T 



(18) 



where P /u " = [</??, </?jj, S' = |[Xi, X°) ^] an d tne CP-odd neutral scalar mass 



matrix is 



By^ + ye 3 ^ + ^ 

"til • ° V\ V\ 



By 



By 



By 

B 3 e 3 ^ + £ 



•V 2 V 2 



-B 3 e 3 



-/^3 
-B 3 e 3 
sl _ £ e a + M 



(19) 



The neutral CP-even scalar sector mass matrix in eq. (18) is given by 



Mg 



(20) 



Byf i +lg 2 v 2 + ye 3 f i +± 



-By - jg 2 z viv 2 



-By, - \g 2 zVl v 2 By% + \g 2 z v 2 - B 3 e 3 ^ + * 



-^3 + I#f^3 

^363 - \g 2 z v 2 v 3 



-^3 + \g\v x v 3 



B 3 e 3 - \g 2 v 2 v 3 — ^ _ R„,„a j_ 1^,2 ^ ta 



where we have defined g| = g 2 + g' 2 . Note that, as a result of CP invariance, the CP-even 
and CP-odd parts of the scalar mass matrices are disjoint and do not mix with each other. 
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The three mass matrices in eqs. flU]) , (|i~9l), and (|20[) are diagonalized by rotation 
matrices which define the eigenvectors 

S + = R S ±S /+ , P = RpoP' , S° = R s oS /0 , (21) 

and the eigenvalues diag(0, m|±, nig±, nig±) = R s ±M|±Rg± for the charged scalars, 
diag(0, nipo, nipo) = R P oMp Rpo for the CP-odd neutral scalars, and diag(m|o, mS^ 2 , m| ) 
= RsoM| RJo for the CP-even neutral scalars. 

The matrices Rs±, Rpo and Rso specify the mixing between the Higgs sector and the 
stau sector. 

If a 3 x 3 matrix M has a zero eigenvalue, then the other two eigenvalues satisfy 

m± = ^TrM ± i ^(TrM) 2 - 4(M n M 22 - M? 2 + M n M 33 - M? 3 + M 22 M 33 - M^) 

(22) 

The CP-odd neutral scalar mass matrix eq. ( |I9"D has a zero determinant, so that its 
eigenvalues m p0 and m 2 p0 (m 2 A and m? fl in the MSSM limit) can be calculated exactly 
with the previous formula. 
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